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ABSTRACT 


The  primary  purpose  of  the  Exhauster  System  of  the  Transonic 
Turbine  Test  Rig  of  the  United  States  Naval  Postgraduate  School, 
Monterey,  California,  is  to  increase  the  range  of  turbine  operating 
pressure  ratios  by  maintaining  a  partial  vacuum  inside  a  hood  that 
encloses  the  turbine.   The  first  part  of  this  study  is  concerned  with 
the  design  and  the  tests  of  the  exhauster  system.   The  experiments 
were  used  to  establish  a  methodlof  predicting  the  operating  character- 
istics of  similar  ejectors  which  is  based  on  the  universal  ejection 
properties  of  turbulent  jets. 

The  second  part  of  this  study  describes  turbine  tests  utilizing 
the  exhauster  system  that  were  carried  out  to  investigate  the  effects 
of  high  values  of  the  isentropic  head  coefficient  and  Reynolds  number 
on  the  turbine  performance.   No  correlation  could  be  established  be- 
tween turbine  Reynolds  number  and  performance  in  the  range  tested. 
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Part  I 
TTR  Exhauster  System  Tests 


1.  Introduction 

The  Transonic  Turbine  Test  Rig  installed  at  the  Propulsion  Labora- 
tory of  the  Department  of  Aeronautics,  U.S.  Naval  Postgraduate  School, 
Monterey,  California,  receives  its  air  supply  from  an  All  is -Chalmers 
VA  312  compressor.   The  maximum  turbine  inlet  total  pressure  attainable 
from  the  compressor  is  about  43  psia.   In  order  to  increase  the  range 
of  possible  operating  pressure  ratios  with  the  existing  air  supply, 
the  turbine  is  enclosed  in  a  hood  which  is  partially  evacuated  by  the 
exhauster  system.   Since  the  turbine  can  only  utilize  a  small  portion 
of  the  mass  flow  from  the  air  supply,  the  remainder  can  be  used  to 
drive  the  exhauster. 

Very  little  design  information  is  available  on  supersonic  ejectors. 
The  first  part  of  this  thesis  is  concerned  with  the  evaluation  of  the 
design  and  performance  of  the  exhauster  system,  and  an  investigation 
of  the  phenomena  associated  with  a  bounded,  supersonic  jet. 

Funds  for  designing  and  building  the  Transonic  Turbine  Test  Rig 
(TTR)  were  furnished  by  the  Bureau  of  Naval  Weapons  (RAPP-14).   Greatly 
appreciated  advice  and  counsel  was  provided  by  Professor  M.H.  Vavra. 
Mr.  L.T.  Clark  of  the  Department  of  Aeronautics  furnished  a  great  deal 
of  advice  and  assistance  in  performing  experiments. 

2.  Installation 

Exhauster  Installation 
The  complete  TTR  and  exhauster  installation  is  covered  in  detail 
by  Eckert.  [s]  To  test  the  exhauster  separately  from  the  turbine,  a 
dummy  exhauster  hood  was  constructed  for  this  purpose  and  secondary 
air  was  taken  in  from  the  atmosphere  through  a  standard  orifice  instal- 
lation.  The  standard  orifice  installation  described  by  Eckert  was 


equipped  with  a  bellmouth  on  the  intake  side  and  a  diffuser  at  the  hood 
entrance.  [4]  Figure  1  shows  a  schematic  of  the  secondary  air  flow  in- 
stallation. 

In  addition  to  the  static  pressure  taps  installed  in  the  exhauster 
and  described  by  Eckert ,  Kiel  probes  were  installed  in  the  nozzle 
plenum  chamber  and  secondary  air  inlet  pipe  (Figure  1).  ["5]  Static 
pressure  tap  station  designations  are  given  in  Figure  1.   All  pressures 
were  measured  in  inches  of  mercury  against  atmospheric  pressure  as  a 
reference,  using  96  inch  manometers  with  .01  inch  graduations,  except 
for  the  driving- nozzle  total  pressure  which  was  measured  on  a  precision 
gauge  graduated  in  increments  of  .05  lbs/square  inch.   Secondary  flow 
rate  differential  pressures  at  the  metering  orifice  were  measured  in 
centimeters  of  mercury  on  a  micro -manometer  graduated  in  increments  of 
.01  cm.   Total  temperature  readings  were  measured  with  standard  I.C. 
thermocouples  using  an  ice  bath  as  a  reference. 
3.   Exhauster-Only  Tests 

General 

Initially  the  nozzle  total  pressure  was  held  constant  while  the 
gate  valve  (Figure  1)  was  used  to  vary  secondary  flow  rates.   Runs 
were  made  with  the  nozzle  total  pressure  set  at  20,  25,  30,  and  35 
psia.   Then  with  the  gate  valve  (Figure  1)  fully  open,  nozzle  total 
pressure  was  varied  from  2  psig  to  the  maximum  obtainable  from  the 
All is -Chalmers  compressor  in  increments  of  2  psig.  Flow  rate  measure- 
ments from  the  standard  orifice,  nozzle  total  pressure  and  temperature, 
and  static  pressure  readings  were  taken  at  each  operating  condition. 
In  order  to  minimize  the  temperature  difference  between  the  primary  and 
secondary  flows,  the  Allis-Chalmers  aftercooler  was  run  at  the  full 
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cold  setting,  giving  a  temperature  difference  of  about  20  degrees  Fahren- 
heit between  the  primary  and  secondary  air.  The  tests  were  then  repeated 
with  one  40  inch  long  section  of  the  mixing  tube  removed. 

Since  the  nozzle  cannot  accommodate  the  full  flow  rate  of  the  Allis- 
Chalmers  compressor,  it  was  necessary  to  throttle  the  compressor  flow, 
discharging  a  portion  of  it  to  the  atmosphere.   The  maximum  total  pres- 
sure obtainable  at  the  driving- nozzle  was  limited  by  approaching  temper- 
ature and  surge  limits  on  the  compressor.   A  six  foot  flexible  pipe 
connection  between  the  driving- nozzle  and  air  supply  (Figure  1)  pro- 
duced an  unexpectedly  large  pressure  loss  (about  8  psig  static  pressure) 
limiting  the  nozzle  total  pressure  to  a  maximum  of  36  psia. 

The  entire  system  was  allowed  to  stabilize  before  taking  readings 
at  each  point  of  operation.   The  system  was  very  stable  except  at  high 
secondary  flow  rates,  where  the  static  pressures  at  stations  one  through 
four  showed  a  large  low  frequency  variation  of  as  much  as  2  inches  of 
mercury.   However,  this  instability  was  not  reflected  in  the  total 
pressure  maintained  inside  the  hood.   The  instability  in  the  mixing 
tube  at  high  flow  rates  is  contrary  to  both  experimental  evidence  and 
theory  given  by  Alexander,  Baron,  and  Comings.  [3]   These  authors  found 
that  the  instability  in  the  mixing  tube  occurred  at  zero  secondary  flow 
rates  and  hypothesized  that  this  was  due  to  the  natural  function  of  the 
jet  in  entraining  air  from  its  surroundings,  requiring  a  reverse  flow 
eddy  in  the  absence  of  any  secondary  flow.   At  zero  secondary  flow 
rates,  the  TTR  exhauster  system  was  very  stable,  even  at  the  maximum 
driving  nozzle  total  pressure  ratio.   It  is  possible  that  hood  leakage 
may  have  furnished  a  slight  amount  of  secondary  flow  to  stabilize  the 
jet  at  zero  secondary  flow  rates,  although  no  evidence  of  any  leakage 
was  detected. 
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The  first  attempt  to  obtain  velocity  profiles  was  made  with  a  W-157 
probe  built  by  the  United  Sensor  Corporation,  East  Hartford,  Conn.,  which 
is  a  wedge  shaped  probe  with  an  included  angle  of  about  34  degrees.   The 
probe  has  a  total  pressure  tap  located  in  a  small  cut-out  section  and 
static  pressure  taps  located  on  the  face  of  the  wedge.   A  mercury  man- 
ometer was  used  to  set  the  probe  to  a  zero  yaw  angle  by  equalizing  the 
static  pressures  from  the  two  face  taps. 

Although  it  was  realized  that  the  detached  shock  Mach  number  for 
the  x^dge  would  be  above  the  Mach  numbers  encountered  in  the  flow,  it 
was  hoped  that  a  calibration  curve  could  be  obtained  over  the  Mach  num- 
ber range  of  interest  (.7  to  1.9).   The  supersonic  wind  tunnel  with 
variable  Mach  numbers  available  at  the  U.S.  Naval  Postgraduate  School 
is  an  AMRAD  Model  W-4  miniature  blowdown  tunnel.   Since  the  physical 
dimensions  of  the  probe  precluded  use  of  the  actual  probe  for  any  cali- 
bration work,  a  full  scale  model  of  the  probe  tip  complete  with  total 
and  static  pressure  taps  was  constructed  so  that  it  could  be  mounted 
in  the  tunnel.   Attempts  to  calibrate  the  probe  on  two  occasions  were 
completely  futile  since  none  of  the  readings  obtained  were  repeatable 
within  a  tolerance  of  i  50  percent. 

Flow  over  the  wedge  was  examined  visually  by  means  of  a  Schlieren 
system  but  failed  to  show  any  unusual  conditions  which  might  have  ac- 
counted for  the  wide  variations  in  dynamic  pressure.   It  is  possible 
that  the  tunnel  calibration  curves  are  inaccurate  or  that  repeatable 
Mach  numbers  cannot  be  obtained  with  the  tunnel.  The  probe  was  tested 
for  leaks,  and  none  were  found  upon  completion  of  testing. 

While  waiting  for  the  construction  of  the  probe  model,  velocity 
traverses  were  taken  with  the  W-157  probe  with  the  concept  that  flow 
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yaw  angles  could  be  checked  even  if  the  probe  could  not  be  calibrated. 
No  evidence  of  any  large  swirl  components  was  found  since  the  yaw 
angles  did  not  exceed  5  degrees  and  were  usually  less  than  1  degree. 
However,  upon  completion  of  the  traverses,  it  was  found  that  the  W-157 
probe  had  developed  a  leak  between  the  total  and  one  static  tap;  there- 
fore, no  conclusions  can  be  based  upon  data  obtained  with  this  probe. 

The  next  attempt  to  take  velocity  traverses  was  made  with  a  pitot- 
static  tube  1/16  inch  in  diameter  mounted  on  a  ^  inch  stainless  steel 
tube.   This  tube  produced  very  erratic  readings  attributable  to  severe 
probe  vibration.   Since  the  probe  did  indicate  a  very  unsymmetrical  flow 
pattern,  this  possibility  was  checked  by  reading  static  pressures  from 
one  orifice  at  a  time  at  each  station.   All  of  the  pressure  taps  pro- 
duced readings  within  .05  inches  of  mercury  of  each  other  showing  no 
indication  of  unsymmetrical  flow. 

Since  the  difficulties  seemed  to  be  due  to  probe  vibration,  a 
total  head  probe  (Figure  2)  was  constructed  so  that  it  could  be  sup- 
ported from  both  ends.   The  idea  of  constructing  a  pitot-static  probe 
was  discarded  initially,  since  the  holes  through  which  the  probe  had 
to  be  inserted  were  too  small  for  a  static  probe  with  the  orifices  lo- 
cated far  enough  back  of  the  tip  to  conform  to  standard  practice  (eight 
to  ten  diameters).   Therefore,  it  was  assumed  that  the  static  pressure 
was  uniform  across  the  section  and  total  head  traverses  were  made  only. 
The  total  head  traverse  at  station  one  indicated  a  lower  total  pressure 
at  the  centerline  than  at  the  off-center  stations.   Since  this  phenomena 
was  not  expected  at  the  time,  a  small  non-standard  static  probe  was 
constructed  anyway  (Figure  2)  to  take  static  traverses.   The  velocity 
profiles  shown  in  Figures  3  to  5  were  obtained  with  these  measurements. 
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No  attempt  was  made  to  calibrate  the  probes  due  to  lack  of  facili- 
ties.  In  addition,  the  probes  could  not  be  inserted  or  withdrawn  while 
the  system  was  in  operation  without  damage  to  the  probe,  so  that  a 
slightly  different  equilibrium  point  was  attained  for  each  station. 

Data  Reduction 

Secondary  flow  rates  were  measured  by  the  orifice  installation  in 
the  secondary  flow  inlet.   Calculation  of  flow  rates  was  accomplished 
by  computer  program  based  on  the  orifice  formulas  used  by  Eckert.  [4] 
Agreement  of  flow  rates  calculated  from  vena  contracta  and  flange  tap 
differential  pressures  were  within  2  percent  for  all  flow  rates  and 
generally  within  1  percent  or  less. 

Nozzle  flow  rates  were  calculated  from  assumed  polytropic  expo- 
nents and  measured  values  of  nozzle  total  temperature  and  pressure. 
Calculations  for  nozzle  flow  rates  in  the  subsonic  region  utilized 
isentropic  flow  and  one- dimensional  normal  shock  relations. 

The  velocity  profile  data  reduction  was  accomplished  by  means  of 
a  computer  program  using  the  Rayleigh  formula,  with  the  known  local 
static  and  total  pressure  behind  the  shock.   Standard  compressible 
flow  relations  were  used  for  subsonic  velocities. 
k.      Exhauster  Performance 

General 

The  classical  method  of  analysis  is  based  on  a  one- dimensional 
momentum  analysis;  however,  if  nothing  is  known  about  the  secondary 
flow  except  its  initial  total  pressure  and  temperature,  the  problem 
cannot  be  solved  since  the  number  of  unknown  quantities  exceeds  the 
number  of  independent  flow  equations  that  can  be  formulated.  There- 
fore, it  became  necessary  to  find  an  additional  relationship  between 
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the  primary  and  secondary  flow  rates. 

Of  all  the  various  approaches  dealing  with  the  problem  of  turbulent 
jets,  the  one  that  appeared  to  have  the  greatest  potential  is  a  hypo- 
thesis set  forth  by  0,  V.  Yakovlevskii.  [14]  His  hypothesis,  based  on 
a  wide  survey  of  experimental  and  theoretical  data, is  that  in  the  dis- 
charge of  a  gas  jet,  the  ejection  properties  are  the  same  as  those  of 
a  geometrically  similar  submerged  jet  of  an  incompressible  fluid  with 
the  same  initial  impulse.   A  submerged  jet  is  defined  as  a  free  jet 
discharging  into  a  fluid  at  rest-   This  hypothesis  leads  to  the  con- 
cept of  the  universality  of  the  ejection  properties  of  a  turbulent  jet 
regardless  of  the  conditions  of  discharge. 

Assuming  that  the  basic  property  of  a  jet  is  its  impulse,  and  that 
the  basic  characteristic  of  such  a  jet  is  the  ejection  or  entrainment 
of  mass  from  the  surrounding  media,  an  arbitrary  jet  and  a  standard, 
isothermal,  submerged  free  jet  discharging  from  geometrically  identical 
nozzles  with  identical  impulses  in  the  initial  cross  section  should 
have  similar  ejection  properties. 

The  basic  part  of  Yakoievskii' s  analysis  is  presented  here  as  a 
basis  for  further  development,  retaining  the  author's  notation: 

Restating  the  basic  hypothesis  in  mathematical  form, 

(1)  oGi'  r  o  Gcj 

dx     dx 

where: 

G  =  gas  discharge,  pVA  in  the  initial  cross  section 

X  =  longitudinal  coordinate  in  the  direction  of  flow 

i  =  index  referring  to  an  arbitrary  jet 

a  =  index  referring  to  a  standard,  submerged,  isothermal  jet 
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Then,  the  relationship  between  the  standard  jet  and  the  arbitrary 
jet  is  expressed  by  the  ratio. 


(2)  -^  ^  V  /  Pcx    : 

^'    V  9; 

o 
where  O^^  is  the  density  of  the  surrounding  medium  and  9  is  the  density 

ratio.  Let: 

(3)  G**  -  — 

Go 

where  Gq  is  the  gas  discharge  at  the  initial  cross  section  and, 


r, 


o 


where  r      is  the   initial   radius  of  an  axisymmetrical   jet.      Eq.    (2)  can 
now  be  written   in  a  non-dimensional   form: 


(5)  dGi   -     rr   de 


dx 


o 


dx 


o 


This  hypothesis  amounts  to  the  assumption  of  the  universality  of  trans- 
verse and  longitudinal  velocity  profiles. 

It  is  shown  by  Yakovlevskii  that  the  increase  in  flow  rate  in  some 
arbitrary  section  of  the  submerged  jet  (Figure  6)  is,  [14] 


(6) 


dG  =  21^  9.-V^/dx 


where : 

0^=  density   in  the  undisturbed  part  of   the  flow 

Voo  =   instantaneous  velocity   in  the  y  direction   in  the  undisturbed 
part  of  the   flow 

Putting   Eq.    (6)    into  a  non-dimensional   form, 
C7)  ^^_      2^     ^    ^,^^ 


d  X  Vo      Uo 


Uvy» 
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where : 

o 

9  =  9**  -  the  ratio  of  the  density  of  the  undisturbed  portion 
"p     of  the  surrounding  flow  to  the  density  of  the  jet 

T  =   radius  at  some  arbitrary  point  of  the  jet 

Ow>  =   instantaneous  centerline  velocity  in  the  x  direction 

<-'oo  =  instantaneous  velocity  in  the  x  direction  in  the  undisturbed 
part  of  the  flow 

•-^o  =  instantaneous  velocity  in  the  x  direction  of  the  jet  discharge 

He  goes  on  to  state  that  the  following  relationship  is  valid  for  this 

case: 

(8)  ^^-         -     ' 


us. 


and  that   for   small   values   of p ,    ^    is   a  constant.      Th 
(9)  ^    Z   -±=r    r  const     -^ 

Then,   using   Eq.    (7)    he  gives   the   formula  for  the   initial   cross- 
section  of   an  axisymmetric    jet   as, 


(10) 


C/^  \/^-'(^"  ^■^.) 


where  ex.,     and  Q,     are  experimentally  determined  constants.   The  follow- 
ing values  are  given  by  Yakovlevskii  for  the  case  of  an  incompressible 
submerged  free  jet:  [l4] 

«i  =  0.11  to  0.20 

^1  =  -4.3 
The  minus  sign  attached  to  o<|  is  due  to  the  Russian  practice  of  basing 
jet  coefficients  on  the  assumption  that  the  jet  is  issuing  from  a  point 
source  located  at  some  distance  in  front  of  the  nozzle,  as  Eq.  (10)  be- 
comes meaningless  near  the  nozzle.   Taking  the  optimum  length  of  the 
mixing  tube  as  8.5  nozzle  diameters,  an  average  of  the  commonly  used 
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value  of   7   to   10   diameters,    and  taking  oC|  as   0.115,   there   is  for  an 
assumed  value  of  V/  P       of  unity, 


9« 


(11)  G°  =  -^ -^'^^Hi^ntfi'^  ''^•^) •  ^-^"^^ 


where: 

4.76  =  driving  nozzle  diameter  (inches) 

8.5   =  optimum  length  of  mixing  tube   (driving  nozzle  diameters), 
taken  as  the  value  for  x 

4.76/2   =   initial   radius  of  the   axisymmetrical    jet,    r^   (inches) 
In  actual ity,\/„£-      varied  from  about  1.2  to  0.96  for  the  present  tests. 

The  optimum  length  chosen  for  the  mixing  tube  was  verified  by  re- 
moving one  UO   inch  section  of  the  mixing  tube.     The  tests  showed  that 
almost      identical  operating  characteristics  were  obtained  for  the  40 
and  80   inch  long  tubes.  Comparing  the  coefficient  of  discharge  of  Eq. 
(11)  with  the  experimental  values  of  the  present  data   (Table  I)   shows 
very  good  agreement  between  the  two  provided  that  a  suitable  choice   is 
made  f or  oC|  •     However,  on  the  basis  of  information  available,   and  due 
to  the  wide  variation  of  o<,  from  author  to  author,   the  arbitrary  choice 
of  o<,i  to   fit  experimental  data  cannot  be  defended.      It   is  noteworthy 
though  that  the  term'G°  is   nearly  a  constant  throughout  the  entire  range 
of  operation  up  to   the   point  of  secondary  flow  choking.      Therefore,    an 
attempt  was  made   to   place  the  flow  relationship  of  Eq.    (10)   on  a  firmer 
basis.      Using  the  standard  definition  of   impulse  function, 

(13)  F  =  PA  ^  p/\V^ 

the  nozzle  impulse  function  was  plotted  versus  \  /  P   where  ^^  is  the 
density  of  the  nozzle  discharge.   Figure  7  shows  that  \/_£_  is  a  linear 
function  of  the  nozzle  impulse  function.   If  the  impulse  function  is 
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plotted  versus  the  secondary  flow  rate  (Figure  8),  a  linear  relation- 
ship is  evident  up  to  the  point  of  secondary  flow  choking,  and  thus  the 
secondary  flow  rate  for  this  particular  system  at  a  given  nozzle  im- 
pulse can  be  determined  from  Figure  8  for  unchoked  secondary  flow.   The 
linear  relation  holds  only  for  the  case  where  the  nozzle  discharge  is 
supersonic.   Since  the  geometry  of  the  nozzle  is  fixed,  there  was  no 
way  to  test  the  generality  of  the  coefficient  of  discharge  determined 
from  Figure  8;  however,  since  the  basic  hypothesis  appears  to  be  true, 
it  should  provide  a  reasonable  first  approximation  for  a  geometrically 
similar  system. 

The  reason  that  the  linear  relationship  between  the  nozzle  im- 
pulse function  and  the  secondary  flow  rate  does  not  hold  for  the  sub- 
sonic case  is  probably  due  to  the  fact  that  one-dimensional  isentropic 
and  normal  shock  relations  were  used  to  calculate  nozzle  exit  condi- 
tions, and  these  are  not  truly  representative  of  the  actual  flow 
phenomena. 

The  next  problem  was  to  find  some  means  of  determining  the  required 
mixing  tube  length.   Since  the  original  hypothesis  was  based  on  the  as- 
sumed universality  of  velocity  profiles,  it  was  assumed  that  a  determi- 
nation of  the  necessary  length  of  the  mixing  tube  could  be  based  on 
incompressible  flow.   However,  even  with  this  simplification,  the  prob- 
lem was  not  tractable  due  to  the  presence  of  a  longitudinal  pressure 
gradient  (Figures  9,  10)  and  the  problems  associated  with  turbulent 
flow.   Therefore,  an  empirical  approximation  was  sought. 

Viktorin  carried  out  a  series  of  experiments  in  1941  using  water 
as  the  driven  and  driving  fluid.  [12]  His  experimentation  was  based  on 
an  earlier  theoretical  analysis  performed  by  Flugel.   Viktorin  found 
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that   the  velocity  distribution   in  the  mixing  chamber  behaved   in  the   same 

1/3 
fashion  as   a  circular  wake  with  the  width  proportional  to   about   x         and 

-2/3 
the  centerline  velocity   proportional   to   about  x  .      The  assumption   is 

made  that   the  mixing  process   reaches   an  optimum  point  when  the  core  of 
the   jet   contacts   the  mixing  tube,   and   that   any  further   increase   in  mix- 
ing tube   length  will   result   in  a  trade  off  between   increased  mass   flow 
rate  and   increased  mixing   and  frictional    losses.      Assuming  that   the 
primary  fluid   emanates   from  a  point   source   located   at    some   point   up- 
stream of  the  nozzle   in  order  to   find   a  constant  of   proportionality, 
it  was  deteirmined  that   the  mixing  process  would  reach  an  optimum  state 
about   40   inches   downstream  of  the   nozzle.      This  was  verified  by  the 
performance   data   (Figures    12,    13  and   Ik)    for  the   system  with  one  40 
inch   section  of   the  mixing  tube  removed. 

A  comparison  of   predicted  centerline  velocities   with  experimental 
data    (Table   II)    shows   that   the   predicted  values   are  within   10%  of  the 
experimental  values   for  stations   two  through  four.      The  discrepancy  at 
station  six  can  be  attributed  to   the  fact   that   at  this   point   the  flow 
is   completely  mixed   and    is   now  characteristic  of  turbulent   flow  in  a 
pipe.      It    is   thought  that   the  flow  at    station  one  differs   from  the   pre- 
dicted value  due  to  off-design  operation  of  the  nozzle,    i.e.,   the 
nozzle  discharge  static   pressure   is   greater  than  the   receiver   static 
pressure,    and  the  nozzle  never   reached   the  design  point  during  the 
experimentation.      Examination  of  the  total   pressure  traverse   shows 
that   the  maximum  total    pressure  was    not  obtained  on  the   axis  of   the 
jet.      This  matches   experimental  data  given  by  Abramovich.    [l]  He  shows 
that  the    jet   actually  has   a  small   subsonic   region  along  the  centerline 
and  that   the  centerline  total   pressure    is   less   than  the  off-axis   total 


20 


pressure  until  the  flow  has  become  sufficiently  adjusted  to  become  iso- 
baric. 

The  analysis   up  to   this   point  has   considered  only  the  case  of  un- 
choked   secondary   flow.      As   is    evident   from  Figure    12,    the   secondary 
flow  rate  becomes   almost   independent   of  the  nozzle  total   pressure   above 
a  pressure   ratio  of   2.3.      A  supersonic   ejector  has   two   regions  of  op- 
eration.     One  region,   the  so-called  mixed   region  occurs   when  the  driving- 
nozzle   total   pressure  ratio    is   sufficiently   low  to  make  the  ejector  flow 
dependent   on  the   back   pressure.      If   the  driving-nozzle   pressure  ratio   is 
increased   sufficiently,   the  flow  will   become   independent   of   the   back 
pressure,    since    it   becomes   supersonic.      As   the  nozzle  pressure  ratio    is 
increased  even  further,    the  secondary  flow  will   eventually  choke   at   the 
point  of  minimum  area. 

The   supersonic  region   is   characterized  by  the  condition  that   the 
secondary   flow  static   pressure   at  the  driving  nozzle  exit   plane   is    less 
than  the  nozzle  exit   static   pressure.      This   point  of  operation  was   not 
reached  with  this   particular  system,    and  the  system  operates    in  the 
mixed   region.      However,    since   secondary  flow  choking  does  occur,    the 
possibility  of   sonic   flow  at   the  minimum  secondary  flow  area  was   checked 
over  a  wide  range  of   nozzle   pressure  ratios,   and    if   was   found  that   the 
secondary   flow  does   not   choke   at   this    point.      Calculation  of   the  poly- 
tropic   exponent   and   the   flow  function    (Equation   14)   between  the   secondary 
flow  inlet   and  the  nozzle  annulus   (Table   III)   shows   that   at   the  higher 
flow  rates,    the  process    is   nearly   isothermal   and   that  the  flow  function 
approaches   the  critical  flow  function    (Equation   17)    indicating  choking 
at   some   point   between  the   inlet   and   the  nozzle  annulus. 

If   the   secondary  flow   is  not  restricted  prior  to   reaching  the 
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nozzle  exit   plane,   flow  characteristics   can  be  found  using  a  method 
developed  by  Fabri.    [t]    A  more  elaborate  method   utilizing  the  method 
of   characteristics   was   presented  by  Addy.    [2l  Both  of  these  methods   are 
concerned  with   flow  in  the  supersonic   region.      For   flow  in  the  mixed 
region,   both  authors   use  a  form  of   the  one-dimensional  momentum  anal- 
ysis.     However,    as  mentioned  before,    solution  of   the  momentum  equation 
depends  on  control  of  the  secondary   flow  rate  or   pressure  ratio. 

The  empirical  method  of  determining  the  ejector  pumping   character- 
istics has   some  serious   drawbacks   as    presented  here.      The  coefficient 
of  ejection  which  is  based  on  empirical   data   is   an  unknown  function  of 
system  geometry,    and  since  the  viscous   effects   have  been   ignored,    it 
does   not   provide  any  means  of   predicting  the  pressure  distribution   in 
the  mixing  tube.      Figures    9   and   10  show  the  mixing  tube  static   pressure 
profiles.      No  theory  could  be  developed  to   produce  the  mixing  tube 
static   pressure  distribution. 

Itomentum  Analysis 

Once  the  values  of   secondary   flow  rate   are  determined   from  Eq. 
(10),    a  one-dimensional  momentum  analysis   can  be   conducted    to  determine 
exhauster  performance.      The  analysis   was  designed  to  parallel  the  pre- 
liminary design  calculations  by  Vavra   in  order  to  determine  the  validity 
of  the  ona— dimensional   approach.     [lOJ    In  order  to   simplify  the  equations, 
polytropic   process   relations   as   well   as   the  non-dimensional   flow  func- 
tion ^  were  used  throughout.      These   were  taken  from  Vavra.    [9] 

For  an  expansion  process, 


<->  ^^VW-  ^  -\m(UlW\ 
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where: 


(15)  Y\'- 


n  =  polytropic  exponent 

?*  =  ratio  of  specific  heats 

W  =  flow  rate  (Ibm/sec) 

T   =  inlet  stagnation  temp.  (OR) 

o 

A  =  flow  area  at  station  of  interest  (in  ) 

p   =  inlet  stagnation  pressure  (psia) 

P  =  static  pressure  at  station  of  interest  (psia) 

R  =  gas  constant  for  air  (ft  Ib/lbm  °R) 

2 
g-  =  universal  gravitational  constant  (Ibm-ft/lb-sec  ) 

For  a  compression  process, 


where  the  subscripts  1  and  2  refer  to  the  inlet  and  station  of  interest 
respectively.   Figure  1  shows  the  exhauster  with  the  station  designa- 
tions used  in  the  analysis. 

For  the  exhauster  operating  conditions  the  flow  through  the  nozzle 
will  be  choked.   For  this  condition,  the  flow  function,^,  becomes 


<->     $. '  i^)\[WT^) 


It  was  assumed  that  the  nozzle  total  pressure  and  temperature  are 
known.   Assuming  a  polytropic  exponent  for  the  flow  in  the  nozzle  to 
the  throat,  ^   can  be  calculated.   Then,  another  polytropic  exponent 
was  assumed  from  the  throat  to  the  nozzle  exit  so  that. 


(18) 


$e  ^$.^ 
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Expressing  the  dimensionless  flow  function  in  terms  of  pressure  ratios, 


"'  *.  ■■  \/iM(^MS^ 


Since  (^    is  known  from  the  choked  flow  condition,    and  P„     and  T„„  are 
known,    it  was   possible  to   solve  for  P,^  by   iteration.      Having  found 

n-> 
(20)  T.N   -  Ttn    ('B^)    -^ 

(21)  ViN  --  \jr^Tc^rvr^rT^\ 


Wn     •■      $e  PtnAnW    -p^ 


(22)  Wn     •■      9eHTNAN\/    -^^ 

Thus  the  complete  operating   conditions  have  been  determined   for 
the  driving  nozzle.      It  was   assumed   that   the  air   leaving  the  nozzle 
would  undergo   a  full   expansion  without  oblique  shocks    in  order  to   con- 
tinue with  the  solution.      If   the  nozzle  total   pressure   is    low  enough, 
the  back  pressure  at   the  exit   plane  will   be  high  enough  to  produce 
shocks    in  the  nozzle  with  attendant   subsonic   flow  at   the  nozzle  exit. 
An   investigation  conducted  on  the  experimental   data  showed  that   the 
flow  from  the  nozzle  would  be   independent  of  the  back  pressure  whenever 
the  nozzle  total   pressure  exceeds    22  psia.      Therefore,   this   particular 
analysis    is  only  valid  for  nozzle  total   pressures   above  22   psia. 

Secondary  air  flow  enters    the  flow  measuring  orifice  with  a  total 
pressure,    PfO'   ^^^  ^  total  temperature,   Trj^Q,    equal   to   atmospheric   con- 
ditions.     Since  the  process    is   essentially  adiabatic,    T-jq  will   remain 
a  constant.      Then  for    the  flow  from   the  hood   to   the  nozzle  exit   plane, 


(23)     W^^/:^  :    $a    Nm  (fer   ■  (fe^o)   "^    ] 
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Because  the  weight  rate  of  flow  of  secondary  air  is  known  for  a  given 
set  of  nozzle  conditions,  a  ratio  of  Pji-'^^TD  ^^^   ^®  assumed,  along  with 
a  value  of  Ptt"   Fo'^  this  assumed  pressure  ratio, 


n-i 


(24)  t,t-Ttd(|^)" 

(25)  VlT    -^29  JCpCTTO-TlT) 

Writing  the  momentum  equation  between  stations   1   and   2,   by  assuming 
uniform  velocity  profiles,    uniform  static   pressure  distributions,    and 
equal    static   pressures   at   the  nozzle  and   secondary   air  discharges, 

(26)(WN2r_WT)V2i_  WnVin-WiVit  ;  A21  ^PiiM-P20.  F-f 

The   friction  force  Ff,   was   expressed   as   a  pressure   loss: 
(27)  F-f  ■-  SLOlii^P  ■-Az\aP 


where, 

(28)  AP  --  -f  i^^N?^ 

Di    2 

L  =    length  of   duct 
D2   =  duct  diameter 
V  =   average  velocity 
p   =   average  density 
f   =   friction   factor 
f    is   a   function  of  Reynolds   number   and   pipe   roughness.      With  an 
approximate  range  of  Reynolds   numbers   from   1.7-2.2   x   10      and  a   surface 
roughness  of    .00015   for   commercial   steel    pipe,    a  roughness   to   diameter 
ratio   of    .0024  was  obtained.      Entering  a  Moody  diagram,    a  value  of  0.017 
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for  f  was  obtained.   Then,  a  aP  was  calculated  and  subtracted  from  P_„ 

IN 

to  obtain  P*. 

Rewriting  equation  (26), 

(29)    /Wfi^WTW2i  -WnVin-WtVit  --Az\(P*-Pz\] 
^  S         /  3  9 

The  total  temperature  after  mixing  was  taken  as  a  weighted  average: 


(30)  Tt2     -    WtTto  +Wn  Ttis 

Wn   +    Wt 

WithWiM^Wr-W   and  writing  the  continuity  equation, 


(31)  Vai  --  WRT2I     r  WR        (Tt2-V2jL\ 

A21  P21         P21A21     \  2Tcp   / 

Solving    (31)    for  P       and  with    (30), 

(32)  Pii    =  WR_   /WtTtd-^WnTtn     _  Vzi^       ] 

A21V21    I  W  ^^^cp   / 

Combining   (30)    and    (32), 

(33)/Wn  -^Wt  V2i)-WnVin-Wi\/it  --Ai\  P*"-  Wr 
\       9  /    g  9  V21 

/  WtTtd   ■»■  WimTtn     •_V2_!____   "\ 

\  W  igJcp      / 

Multiplying  (33)  by  ¥2;)^  and  expanding, 


(34)    V2|^ 


Let: 

(35) 

(36) 


(Wn-^Wt   _  WR         "i    -  V21  fWN  ViN   ^WtVit 
\       g  2gTcp       /  ^  9 

■^  A2>  P*\     t   WRT    -.0 

B 1  =  W  h  -_B ) 

9     \       2gTcp      / 
82    ^      -  /  Wisi  ViN  -^  Wt  ViT     i-PiAz\ 
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(37)  B3    =  WRTt2 

Solving  the  quadratic    (34)    for  velocity, 


(38)  V21    --    -.^2.    +\//Bxf -El2 

2  61     ~]j{Tq7}      Bi 
This  equation  yields   both  a   subsonic  and   a  supersonic   solution; 
however,   only  the  subsonic   solution  was   pursued   since  the   supersonic 
solution  did  not   provide  sufficient   pressure   recovery  under   the  opera- 
ting conditions   investigated.     With  the  known  quantities, 

(39)  Tzi   --  Tt2  -  Vzi^ 

ZcjJcp 

±zL 

(40)  Pt2    ■-    PzJTtZ  ^  ^ 


(^) 


For  the  subsonic   solution,   the   flow  undergoes    a  compression  pro- 
cess   in  the  diffuser: 


(41)  w\/t2i\/r"  --  (f) .-  \/2?-rTT2  ip^]^^'-(p^y^] 

This  equation  was  solved  by  iteration  to  find  P3.   If  P3  does  not  match 
^AMB  (atmospheric  pressure)  a  new  value  for  P-rj^  was  assumed  and  the 
process  repeated  until  P3  matched  P^^j^- 
Having  obtained  P3, 

(^2)  T3:    Tz>/P3  \"^ 

(^3)  VS  ^   \/2^1cp(TTl-Tl\) 


(44) 


PT3     --     P3/Tt2_\ 
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Test  Results 

The  equations  just  described  were  used  to  predict  exhauster  per- 
formance utilizing  a  computer  program.   Since  the  primary  purpose  of 
the  exhauster  is  to  maintain  a  partial  vacuum  inside  the  hood,  this  was 
the  major  point  of  interest.   A  comparison  of  experimental  and  theoret- 
ical data  is  given  in  Figure  11.   The  discrepancy  is  due  to  the  fact 
that  the  one-dimensional  analysis  ignores  conditions  from  the  inlet  to 
the  hood,  assuming  that  the  polytropic  exponents  are  constant.  This 
is  not  the  case  since  secondary  flow  conditions  are  dependent  on  the 
varying  polytropic  process  from  the  inlet  to  the  hood.   Table  III  shows 
the  variation  in  polytropic  exponents. 

For  the  case  where  the  secondary  flow  into  the  hood  is  a  constant 
(i.e.,  turbine  operation),  the  one-dimensional  momentum  analysis  ap- 
pears to  provide  a  good  approximation  (Figure  11).   It  is  to  be  noted 
that  the  experimentation  was  conducted  with  a  nozzle  total  temperature 
about  20Of  above  the  secondary  flow  total  temperature.   Increasing  the 
nozzle  total  temperature  will  result  in  a  very  slight  decrease  in  ex- 
hauster pumping  capability,  and  for  the  allowable  range  of  turbine 
operating  temperatures,  exhauster  flow  is  virtually  independent  of 
total  temperature. 
5.   Exhauster  Design  Evaluation 

General 

The  most    recent  work  on  the  design  of  ejectors   available  was    pub- 
lished  by  Engle   in   1963.    [6]  He   presents   a   set  of  equations   to   be  used 
to   optimize  the   design  geometry  based  on  a  one— dimensional   analysis   of 
incompressible  flow.      The  present  exhauster  system  meets   the  criteria 
set    forth  by  Engle   for  driving-nozzle   placement   and   secondary  flow 
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entry.   For  optimum  operation,  the  two  streams  should  be  parallel  upon 
mixing  tube  entry  to  minimize  losses.   Secondary  flow  static  pressure 
profiles  (Figure  14)  show  a  continuous  decrease  in  static  pressure  up 
to  the  nozzle  exit  plane  except  at  very  high  nozzle  pressure  ratios. 
The  increase  in  static  pressure  at  station  28  at  high  nozzle  pressure 
ratios  is  probably  due  to  separation. 

At  high  secondary  flow  rates,  the  mixing  tube  static  pressures 
show  a  fluctuation  of  1.0  to  2.0  inches  of  mercury.   However,  the  flow 
rate  and  hood  total  pressure  remained  steady  in  this  regieme  of  opera- 
tion.  Since  the  secondary  flow  is  choked  at  these  pressure  ratios,  it 
was  thought  that  the  secondary  flow  rate  might  not  be  sufficient  to 
stabilize  the  jet.   This  was  not  the  case  though  for  this  system.   The 
exhauster  operation  was  very  stable  at  high  nozzle  pressure  ratios  with 
no  secondary  flow.   The  instability  only  appears  at  high  pressure  ratios, 
without  being  felt  inside  the  hood.   The  probable  cause  of  the  instabil- 
ity is  the  formation  of  a  recirculation  eddy  due  to  separation  effects 
along  the  mixing  tube  since  the  secondary  flow  is  under  an  adverse 
pressure  gradient. 

Engle  flOj  also  presents  a  method  for  determining  the  required  area 
of  the  mixing  tube.   This  is  based  on  continuity,  and  is, 


(^5)  A  --  \a/ 

Vai  ViN  9 

Substituting  experimentally  determined  values    into  this   equation  for   a 


nozzle  pressure   ratio  of   2.22, 


(46)  A   -■      (  8.0I1   ^  3.gi       --.68Gf+* 
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where  ^  is  an  average  density  based  on  the  average  velocity  computed  at 
station  five.   The  actual  area  of  the  mixing  tube  is  .306  ft^.   Since 
Engle  states  that  in  practice  the  mixing  tube  area  should  be  some  30  to 
50  percent  larger  than  the  calculated  value,  the  mixing  tube  area  of 
this  system  is  about  30%  of  the  desired  optimum. 

Optimum  mixing  occurs  according  t6' Engle  when  the  ratio  between 
the  velocity  in  the  core  and  in  the  annulus  is  between  .7  and  .8  and 
that  the  length  required  is  a  function  of  the  nozzle  to  secondary  flow 
area.  ^6^  In  general,  the  optimum  mixing  length  is  between  7  and  10 
diameters  of  the  mixing  tube.   Removal  of  one  section  of  the  mixing 
tube  left  the  exhauster  operation  virtually  unchanged  except  at  the 
higher  pressure  ratios.   The  velocity  profiles  taken  at  the  diffuser 
entrance  (Figure  5)  for  the  short  mixing  tube  indicate  that  the  ratio 
between  the  core  and  annulus  velocity  increases  with  increasing  nozzle 
pressure  ratio,  and  falls  short  of  the  . 7  to  .8  desired  velocity  ratio 
which  accounts  for  the  reduction  in  pumping  capability  at  the  higher 
pressure  ratios. 

Optimum  mixing  tube  length  for  this  particular  exhauster  falls 
between  the  short  and  long  tubes.   However,  as  indicated  from  the  ex- 
hauster performance,  a  wide  variation  in  mixing  tube  length  has  little 
effect  on  system  performance,  since  varying  the  length  to  diameter  ratio 
from  5.46  to  10.9  had  only  a  small  effect  on  the  operating  characteris- 
tics. 

It  was  not  possible  to  duplicate  turbine  operating  conditions  as 
far  as  secondary  flow  rates  and  nozzle  total  pressure  are  concerned. 
However,  extrapolating  from  experimental  data,  it  should  be  possible 
to  maintain  a  pressure  ratio  of  about  5  to  1  for  the  converging -diverging 
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stator  and  a  pressure  ratio  of  about  6.3  to  1  for  the  converging  turbine 
rotor  stator  for  the  transonic  turbine.   Exhauster  operating  conditions 
can  be  determined  from  Figures  11,  12,  and  13. 

Discussion 

The  basic  hypothesis  set  forth  by  Yakovlevskii  appears  to  be  valid 
since  the  ejection  property  of  a  given  nozzle  is  a  function  of  its  im- 
pulse, and  empirical  relations  found  for  incompressible  flow  can  be 
used  to  find  a  first  approximation  for  the  case  of  compressible  flow. 
If  the  secondary  flow  rate  or  pressure  is  known,  a  one-dimensional 
momentum  analysis  can  be  used  to  determine  operating  conditions  pro- 
vided that  the  secondary  flow  is  not  choked.   If  the  secondary  flow 
is  choked  at  the  mixing  chamber  inlet,  the  method  of  Fabri  can  be  used 
to  predict  operating  characteristics.   The  method  based  on  the  hypothe- 
sis of  Yakovlevskii  has  the  advantage  of  producing  a  good  approximation 
of  the  pumping  characteristics  with  relative  ease.   However,  the  gen- 
erality of  the  relationship  between  the  secondary  flow  and  nozzle  im- 
pulse function  was  not  tested  for  other  ejector  configurations. 

Determination  of  the  required  length  of  mixing  tube  for  a  given 
ejector  is  still  largely  a  matter  of  experience^  however,  it  was  shown 
that  a  wide  variation  in  mixing  length  had  very  little  effect  on  system 
operation.   Mixing  tube  length  to  diameter  ratios  of  from  7  to  10 
should  produce  optimum  operation.   Calculation  of  the  required  mixing 
tube  diameter  based  on  the  equation  given  by  Engle  indicates  that  the 
mixing  tube  diameter  falls  short  of  the  desired  optimum. 

Although  unstable  conditions  did  appear  in  the  mixing  tube  at  high 
flow  rates,  the  instability  was  not  reflected  in  the  hood  total  pres- 
sure.  This  fact  coupled  with  the  extrapolated  performance  figures 
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indicates  that  the  exhauster  will  provide  turbine  operating  pressure 
ratios  in  the  design  range. 
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Part  II 
Tests  of  a  Reaction  Turbine 
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1.  Introduction 

The  second  portion  of  the  thesis  describes  a  turbine  test  program 
carried  out  with  the  exhauster  system  to  investigate  the  effects  of 
high  values  of  the  isentropic  head  coefficient  on  turbine  performance. 
The  turbine  tests  were  used  also  to  correlate  turbine  performance  with 
Reynolds  number. 

2.  Turbine  Installation 

The  turbine  installation  is  identical  to  that  described  by  Eckert 
except  as  noted  here.  [Sj  The  turbine  is  an  ARES  MOD  II  instead  of  the 
Transonic  Turbine.   A  cross  section  of  the  turbine  is  shown  in  Figure 
15.  Figures  16  to  18  show  details  of  the  rotor  and  stator  construction. 
The  instrumentation  was  modified  by  placing  additional  pressure  taps 
in  the  shroud  as  shown  in  Figure  15.  .  It  is  to  be  noted  that  the  pres- 
sure tap  locations  have  been  projected  into  the  one  dimentional  view 
while  in  actuality  they  are  spaced  around  the  periphery  of  the  shroud. 
For  this  particular  test  program,  the  axial  clearance  was  set  at  0.41 
inches  while  the  radial  clearance  was  set  at  0.033  inches. 

The  cover  plate  flexure  instrumentation  mentioned  by  Eckert  was 
installed  to  permit  measurement  of  the  torque  and  the  axial  force  on 
the  cover  plate.  [5]  - 

3.  Turbine  Test  Program 

Flow  Nozzle  Calibration 
The  turbine  flow  nozzle  installation  and  calibration  techniques 
used  are  covered  by  Eckert.  [4]  Early  calibration  runs  conducted  by 
Eckert  and  Mr.  L.  T.  Clark  indicated  that  at  low  flow  rates,  the  nozzle 
coefficient  seemed  to  be  a  function  of  supply  total  pressure-   Since 
this  violates  the  principle  of  similarity,  additional  calibration  runs 
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were  made  at   supply  total   pressures   of   20,    22,    and   24  psia.      Data  was 
reduced  using  computer  program  FLOCAL  with  the  exception  that   the  ex- 
pansion  factor  Y-i    for  the   nozzle  was   changed  to   that  given  by  the  ASME 
Power  Test   Codes.    [SJ    Figure   19   shows   the  nozzle  coefficient   as   a  func- 
tion of  Reynolds  number.      The  nozzle  coefficient  was   based  on  vena 
contracta  tap  data  since  these  have   a   lower  tolerance.      Examination  of 
earlier  calibration   runs  reveals    a   fairly   large   scatter   (as  much  as   ten 
percent)   between  vena  contracta  and   flange  tap  data  at    low  flow  rates. 
For  the  calibration   run  shown  in  Figure   19,   the  difference   between  vena 
contracta  and  flange  tap  calculated   flow  rates    is    less   than  one  per- 
cent.     Therefore,   a  nozzle  coefficient  based  on  Figure   19  was  used.      An 
analytical   expression  for  the  flow  nozzle  coefficient   as   a   function  of 
Reynolds    number  was   found  by  using  the  method  of    least   squares   to  ob- 
tain a  polynomial    approximation.      This  was    accomplished  by   a  computer 
program  yielding  an  eighth  order  polynomial   for  the  best   fit.      The   poly- 
nomial   is, 

C   =  -69.4248    +    65.42325X   -   25.4166x2    +   5.38697x3 
(47) 

-0.67816x^   +  0.051422x5   -  0.0022472x6   +  0. 000049 356x'' 

-0.00000035959X® 

where  x  =  Reynolds   number  multiplied  by    10"^.      The  maximum  deviation 
in  the  flow  range  of   interest    (1.6   to  4.5    Ibm/sec)    occurs   at   a  Reynolds 
number  of   6  x   10^   and    is    .003.      The   polynomial    is  only  valid   for  Rey- 
nolds  numbers   in  the   range  of  4  to    12   times    10    . 

Turbine  Flow  Rate  Determination 
The  turbine  flow  rate  is   equal   to  the  total   rate  measured  by  the 
nozzle   less   the  plenum   labyrinth   leak  rates.      Total   nozzle  flow  was 
determined   by  means   of   the  ASME   standard  equations   using   an   initial 
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nozzle  coefficient  of  1.03.  ^15]  Vena  contracta  tap  data  was  used 
throughout,  with  differential  pressure  measured  on  a  96  inch  water 
manometer,  and  flow  nozzle  static  pressure  measured  on  a  96  inch  mer- 
cury manometer.  Then,  after  calculating  the  Reynolds  number,  the  flow 
nozzle  coefficient  was  corrected  by  application  of  the  eighth  order 
polynomial  described  above  and  a  new  flow  rate  determined. 

Calculation  of  the  plenum  labyrinth  leakage  flow  rates  utilized 
the  expressions  given  by  Eckert  which  were  determined  from  experimental 
data.  [4]  . 

Turbine  Tests 

The  turbine  was  rotated  on  14  November  1966  with  the  bottom  portion 
of  the  hood   installed  to   check  the   instrumentation.      The  shaft  seal 
labyrinth  was  bored  out  to   1.258   inches  diameter  to  accommodate  a  new  " 
quill  shaft  without  recutting  the  chambers   for  a  seal  tooth  clearance 
of  0.040   inches.      The  turbine  was   run  up  to   about   IJ, 000  RPM.      Upon 
shutting  down,    it  was   noticed  that,  the  quill   shaft  had  rubbed  on  the 
labyrinth,   scoring  the  shaft.      The  shaft  was  buffed  and  reinstalled 
with  the   labyrinth  realigned.      When  the  turbine  was   rotated   again  on 
16  November,    shaft   rubbing  occurred   at. about    7,000  RPM.      The  shaft  was 
measured   and   found  to  be  0.005   inches  out  of   round.      It  was   then  turned 
down  0.005    inches   and  the  labyrinth  enlarged  so  that   after  alignment 
thfe  clearance  was   0.030   inches   radially  at   the  dynamometer  end   and  0.015 
Inches   radially  at  the  hood  end.      This  prevented  further  rubbing  when 
the  turbine  was  rotated  on  18  November  1966  with  the  full  hood  and  ex- 
hauster  installed.      The  present   labyrinth   is  only  a  temporary  device 
and  will   be  replaced.      Since  the  only  result   of  the  increased   leakage 
Into  .the  hood   is   the  degradation  of  the   system's   capability  to  maintain 
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a  partial  vacuum  inside  the  hood,  the  existing  labyrinth  is  suitable 
for  continued  testing  until  a  replacement  is  obtained.   The  purpose  of 
the  first  run  was  to  Investigate  the  performance  limits  of  the  system. 
With  the  ARES  MOD  II  turbine  installed,  a  pressure  ratio  of  1.79  was 
attained.   Since  the  turbine  efficiency  is  much  greater  than  that  of 
the  exhauster,  a  reduction  in  turbine  flow  rate  with  the  corresponding 
increase  in  exhauster  flow  rate  resulted  in  a  decrease  in  pressure  ratio. 
The  exhauster  is  capable  of  maintaining  a  partial  vacuum  of  one-half  an 
atmosphere  in  the  hood  with  a  turbine  flow  rate  of  3.88  Ibm/sec.   Since 
this  particular  turbine  has  been  the  subject  of  an  extensive  test  pro- 
gram without  the  hood  installed  reaching  a  maximum  pressure  ratio  of 
1.55,  it  was  decided  to  run  at  pressure  ratios  of  1.3,  1.4,  1.5,  1.55, 
and  1.65  with  values  of  the  isentropic  head  coefficient  varying  between 
the  design  value  of  2.609and  4.5. 

Turbine  testing  was  started  on  22  November  1966.   During  the  run, 
the  RPM  readout  was  checked  by  means  of  a  strobotac  and  it  was  found 
that  the  RPM  read  out  was  faulty.   The  system  was  checked  thoroughly, 
and  the  difficulty  finally  corrected  by  replacing  the  leads  between  the 
counter  and  the  flux  cutter  at  the  turbine.   The  final  test  was  con- 
ducted on  23  November  1966,  when  testing  was  halted  due  to  time  con- 
siderations.  Under  the  present  test  rig  set-up,  the  flux  cutter  and 
pickup  is  located  inside  the  hood.   It  is  felt  that  it  would  be  advan- 
tageous to  move  the  unit  outside  the  hood  where  it  would  be  readily 
accessible  for  adjustment. 

Data  Reduction  and  Results 

The  turbine  flow  rate  was  determined  as  described  in  the  previous 
section.   Total  inlet  temperature  was  determined  from  a  probe  installed 
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in  the  6  inch  pipe  of  the  stator  assembly.   Inlet  total  pressure  was 
taken  from  the  average  of  the  six  Kiel  probes  (Figure  15)  placed 
radially  at  the  stator  inlet.   The  hood  pressure  taken  from  the  average 
of  the  eight  hood  static  pressure  taps  was  used  for  the  turbine  discharge 
static  pressure.   The  turbine  speed  and  dynamometer  moment  complete  the 
quantities  required  to  determine  turbine  performance.   Data  was  reduced 
using  the  method  of  Vavra  by  a  computer  program.  The  method  of  analy- 
sis is  being  published  in  a  Thesis  by  Lt.  M.  W.  Wallace,  U.S.  Naval 
Postgraduate  School,  lyfonterey,  California.  [l3j 

The  following  referred  values  were  obtained  from  the  performance 
analysis ; 

for  y  =cp/cv  =  1.4: 

(49)  ©  =   Tto  /5I8.H 

(50)  <f  =   Pro/ I M.  7 

(51)  M  rci      =r^  =  referred  moment  (Ft-lb) 

(52)  N  vtf   "JNr  =  referred  turbine  speed  (RPM) 

S 

(53)  Wvtf       ^WV*©"  ~  referred  flow  rate    (Ibm/sec) 

~i 

(54)  upv-tf     =HP       -  referred  horsepower   (HP) 

\7©T 

(55)  K  I  s         =  ^ ^ ; s     =  head  coefficient 

U.V2 

The  turbine  efficiency  is  a  total  to  static  efficiency,  defined  by 


(56)  }^     =±LP_ 

HPth 


^  V© 

where  HP^r^   is   the  theoretical  horsepower  based  on  an   isentropic   expan- 
sion from  the   inlet  total  temperature   to   the  static   pressure   in  the 
hood ,   or , 
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r  -^ 

(57)  HPtv,         =/W\/©'\       I     -/P2    ^^ 


70- '3987 


where 

p2  =  discharge  pressure    (Ib/ft^) 

Pro  =   inlet   total   pressure    (Ib/ft^) 
Determination  of  turbine  Reynolds   numbers    is   rather  arbitrary  as 
to   choice  of  characteristic    lengths    and  velocities.      After   a  perusal 
of   the   literature  on  the  subject  of  Reynolds   number   effects    in  turbo- 
machines,   two   sets  of  quantities   were  chosen,   one  based  on  rotor  tip 
speed,    and  one  based  on  the  flow  rate  and  mean  radius.     They  are  de- 
fined  as, 

(58)  RE  -  W_  .     ' 

oAAYvy 

and 

(59)  RT  ^  ZIY  RPMp  R 

120  ju 
where 

W  =  weight   rate  of  flow   (slugs/sec) 

2 
^AA  =  viscosity  at  turbine  inlet  total   conditions    (lb. sec/ft   ) 

(>    =  density   at   inlet  total   conditions    (slugs/ft^) 

R  =  rotor  radius    (inches) 

^y^  =  mean  rotor  radius    (ft) 

The  referred  values  of  horsepower,   moment,    and   flow  rates   are   shown 

in  Figures   20  through  22.      Turbine  efficiency  plotted   as   a  function  of 

the   isentropic  head  coefficient    is   shown   in  Figures    23  and  24.      Although 

great   care  was   used   in  attempting  to   keep  the  pressure   ratio   at    the  set 

t 
value,   the  wide   number  of    system  variables   caused  minor  deviations   from 

the  ratio.      Table   IV  gives  the  actual   pressure  ratio   at  each  data  point 
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and  the  maximum  deviation. 
4.   Discussion  of  Results 

General 

The  exhauster  system  functioned  well,  maintaining  a  hood  pressure 
of  about  one  half  of  the  atmospheric  pressure.   At  a  pressure  ratio  of 
1.65  with  a  turbine  flow  rate  of  3.88  Ibm/sec  a  bothersome  low  fre- 
quency instability  was  noted  in  the  hood  pressure  which  was  reflected 
in  turbine  operation.   This  fluctuation  resulted  in  a  low  frequency 
flow  rate  fluctuation  of  as  much  as  two  inches  of  water  in  differential 
pressure  across  the  flow  nozzle  and  an  RPM  variation  on  the  order  of 
150.  All  of  the  turbine  pressure  readings  were  made  on  one  manometer 
bank  by  using  a  Polaroid  camera  to  photograph  the  board,  and  since  the 
flow  rate  differential  pressure,  turbine  RPM  and  dynamometer  reading 
were  taken  as  quickly  as  possible,  the  fluctuation  did  not  seem  to  in- 
duce undue  scatter  in  the  data. 

The  torque  capsules  are  calibrated  prior  to  each  run  by  applying 
known  v/eights  to  a  lever  arm.   For  dynamometer  calibration,  the  torque 
capsule  must  be  removed  from  the  unit  and  placed  in  a  special  stand. 
The  present  calibration  technique  involves  jiggling  or  tapping  the  cap- 
sule as  the  weights  are  varied  through  several  cycles  until  the  read- 
out produces  consistent  values.   For  the  last  run,  the  stator  axial 
force  capsule  could  not  be  calibrated  at  all,  but  since  data  reduction 
by  the  method  used  does  not  require  this  value,  no  attempt  was  made  to 
correct  this  difficulty.   It  is  felt  that  the  present  method  of  cali- 
brating the  torque  capsules  mechanically  is  unsatisfactory  both  from 
the  standpoint  of  accuracy  and  time.  Calibration  of  the  torque  cap- 
sules required  as  much  as  three  hours,  severely  limiting  the  amount  of 
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experimentation  that   could  be   accomplished  during  the  normal   working  day. 
Provided  that   the  capsules   are  not  disturbed,    electrical   calibration 
would  be  both  rapid   and   accurate.      The  present   axial   force  measurements 
when  used   in  a  momentum  analysis  do   not   produce  results   consistent  with 
those  obtained  from  continuity,    and  the  difficulty  almost   certainly   lies 
with  the  axial   force  measurements. 

The  turbine   is  designed   for   an   isentropic  head  coefficient  of  2.609 
at   a   pressure   ratio  of    1.5.      Neither    increasing  the   pressure   ratio   nor 
increasing  the   isentropic  head  coefficient  did  result    in  an   increase   in 
efficiency    (Figures   2  3  and   24).      At   pressure  ratios  of   1.3  and   1.65   a 
dip   in  efficiency  was   noted  with   increasing  head  coefficients   that    is 
not  present   at  the  other  pressure  ratios.      The  reason  for  this  dip   is 
unknown,   and   it  does   not   match  the  theoretical   data  given  by   Vavra.^lll 
It    is   to   be   noted  that   this   particular  turbine   is    capable  of  operating 
under  off-design  conditions  with  very   little  reduction  in  efficiency. 

Although  a  great  deal   of  research  has  been  devoted  to   investigating 
Reynolds   number  effects    in  turbomachines ,   the  results   are   still    incon- 
clusive.     Figure  25  shows   the   Reynolds   number  based  on  flow  rate  plotted 
versus  turbine  efficiency  for  three  runs.     Run  50   was   conducted  by  the 
author  using   the  full  hood  and  exhauster   installation.      Runs   40  and  45 
were  conducted  by  Professor  M.H.    Vavra  with  the  turbine  discharging  to 
the   atmosphere   at   pressure  ratios  of   1.3,    1.4,    and   1.5.      Run  45   used 
the  same  axial   clearance  of  0.41    inches    and   tip  clearance  of   0.0  33 
inches   that   were  used  on  Run  50,   whereas   for  Run  40,    at   an   axial   clear- 
ance of  0.410    inches   the   tip  clearance  was   set    to   0.015   inches.      For  the 
plotting  of   the  data,    the  efficiency   at  each  pressure  ratio   was   averaged 
and  plotted   as   a  single   point   since   the  Reynolds   number  variation  at   a 
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given  pressure  ratio  is  insignificant.   Reynolds  number  based  on  rotor 
tip  speed  is  shown  plotted  versus  turbine  efficiency  in  Figure  26. 

NASA  investigators  on  a  number  of  turbines  using  the  flow  rate 
definition  of  Reynolds  number  found  that  turbines  operating  with  Rey- 
nolds numbers  above  2  x  10^  showed  no  variation  in  loss  parameters  with 
Reynolds  numbers.  fS]  They  also  found  that  different  turbines  operating 
in  the  same  range  of  Reynolds  numbers  had  different  loss  parameters. 
This  seems  to  point  out  that  either  the  losses  in  turbines  operated  above 
a  Reynolds  number  of  2  x  10  are  not  a  function  of  Reynolds  number  or 
that  the  machine  characteristics  used  to  determine  the  Reynolds  number 
are  not  correct. 

Examination  of  Figure  26  which  uses  the  Reynolds  number  based  on 
rotor  tip  speed  fails  to  show  any  trend,  and  again  it  appears  that  the 
efficiency  is  not  a  function  of  the  Reynolds  number  within  the  range  of 
investigation. 

Losses  in  turbomachines  are  due  to  viscous  effects  and  since  Rey- 
nolds number  is  a  measure  of  these  effects,  it  vrould  seem  only  natural 
that  some  correlation  should  exist  between  Reynolds  number  and  turbine 
efficiency.   The  primary  difficulty  lies  in  the  elusiveness  of  the 
actual  loss  mechanisms  associated  with  the  flow  in  turbomachinery,  and 
therefore,  the  inability  of  the  investigator  to  be  capable  of  relating 
the  loss  parameter  to  Reynolds  number.   A  more  detailed  study  would  in- 
volve various  definitions  of  characteristics  lengths  and  velocities  over 
a  wide  range  of  Reynolds  numbers.   This  is  a  difficult  task  experimentally 
when  dealing  with  a  compressible  fluid  due  to  the  wide  range  of  tempera- 
tures and  pressures  required  to  obtain  a  significant  variation  in  Rey- 
nolds number. 
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5.   Conclusions  and  Recommendations 

The  ARES  MOD  II  turbine  appears  to  reach  its  maximum  efficiency  at 
the  design  pressure  ratio  and  the  design  value  of  the  isentropic  head 
coefficient.   However,  it  is  capable  of  operation  over  a  wide  range  of 
pressure  ratios  and  isentropic  head  coefficients  with  only  a  slight  re- 
duction in  efficiency.   For  the  operating  ranges  tested,  turbine  per- 
formance does  not  seem  to  be  related  to  Reynolds  number. 

The  weakest  point  of  the  present  turbine  test  rig  is  the  torque 
capsule  calibration.   A  suitable  means  of  electrical  calibration  should 
provide  both  speed  and  accuracy  in  calibration.   The  axial  force 
measurements,  both  on  the  stator  and  coverplate,  seem  to  be  in  error, 
due  to  some  unknown  cause.   Therefore,  it  is  recommended  that  the  pre- 
sent instrumentation  be  improved,  possibly  by  installing  some  type  of 
hydraulic  measuring  device. 
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Table   I 


C?       as   a   function  of    Pt  Nome  /  P/>mb 

"F"  Ptnoj:7up/Pamb 

1.481  i.ys"*" 

1.482  1.88"^ 
1.482  2.10'^ 
1.478  2.16'^ 
1.468  2.36"^ 

1.482  i-ys"^-^ 

1.480  2.10'^"^ 

1.480  2.22'^"^ 

1.470  2.50'^'^ 

1.421  1.14'''^"^ 

1.468  1.29'^"^'*' 

1.482  1.41"^'^'^ 


+      Long  Mixing  Tube 

++     Short  Mixing  Tube 

+++  Subsonic   Nozzle  Operating  Region 
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Table  II 
Centerline  Velocities,  ft/sec,  Ptn/Pamb  =  2.22 

Station  Theoretical  Experimental 

1  1400  1100+ 

2  .   "  1270  1260 

3  920  960 

4  7  35  760 

6  570  650++ 

+  Station  1  is  the  station  where  the  centerline  total  pressure 
is  much  lower  than  the  off-center  total  pressure,  and  the  experi- 
mental value  is  in  question. 

++   Station  6  is  located  sufficiently  far  downstream  of  the 
point  where  the  mixing  process  is  complete  and  the  flow  is  that  of 
turbulent  flow  in  a  pipe  and  is  not  expected  to  c6nform  to  the 
power  law  assumed  for  the  centerline  velocity. 
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Table   III 
<J  and    <J> critical   ^®   ^  Function  of  Secondary  Mass   Flow  Rate 
Process  from  Secondary  Flow  Inlet  to  Annulus 


* 

^critical 

W(lbm/sec) 

n 

2614 

.2667 

3.43 

1.056 

2755 

.2759 

3.62 

1.060 

2876 

.2891 

3.785 

1.066 

2986 

.3119 

3.93 

1.0770 

3046 

.3348 

4.00 

1.089 

Process   from  Hood  to  Annulus 


$ 


.3509 
.3952 
.4475 
,5150 
.5616 


ritical 

W(lbm/sec) 

Y\ 

5950 

3.43 

1.295 

5874 

3.62 

1.287 

6015 

3.78 

1.302 

60  34 

3.93 

1.304 

6125 

4.00 

1.3140 
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